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Abstract 

As shown in previous work, there is a well-defined nonperturbative gravitational 
path integral including an explicit sum over topologies in the setting of Causal Dy- 
namical Triangulations in two dimensions. In this paper we derive a complete ana- 
lytical solution of the quantum continuum dynamics of this model, obtained uniquely 
by means of a double-scaling limit. We show that the presence of infinitesimal worm- 
holes leads to a decrease in the effective cosmological constant, reminiscent of the 
suppression mechanism considered by Coleman and others in the four-dimensional 
Euclidean path integral. Remarkably, in the continuum limit we obtain a finite space- 
time density of microscopic wormholes without assuming fundamental discreteness. 
This shows that one can in principle make sense of a gravitational path integral which 
includes a sum over topologies, provided suitable causality restrictions are imposed 
on the path integral histories. 
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1 Introduction 



Despite recent progress PI Ej, little is known about the ultimate configuration space of 
quantum gravity on which its nonperturbative dynamics takes place. This makes it diffi- 
cult to decide which (auxiliary) configuration space to choose as starting point for a quan- 
tization. In the context of a path integral quantization of gravity, the relevant question is 
which class of geometries one should be integrating over in the first place. Setting aside the 
formidable difficulties in "doing the integral" , there is a subtle balance between including 
too many geometries - such that the integral will simply fail to exist (nonperturbatively) 
in any meaningful way, even after renormalization - and including too few geometries, with 
the danger of not capturing a physically relevant part of the configuration space. 

A time-honoured part of this discussion is the question of whether a sum over different 
spacetime topologies should be included in the gravitational path integral. The absence to 
date of a viable theory of quantum gravity in four dimensions has not hindered speculation 
on the potential physical significance of processes involving topology change (for reviews, 
see OH]). Because such processes necessarily violate causality, they are usually considered 
in a Euclidean setting where the issue does not arise. Even if one believes that Euclidean 
quantum gravity without a sum over topologies exists nonperturbatively as a fundamental 
theory of nature - something for which there is currently little evidence -, insisting on 
including such a sum makes it all the more difficult to perform the path integral and 
extract physical information from it. This happens because the number of ways in which 
one can cut and reglue a manifold to obtain manifolds with a different topology is very 
large and leads to uncontrollable divergences in the path integral.^ 

Formal path integrals involving nontrivial topologies are either semiclassical and as- 
sume that a (usually small and hand-picked) class of such configurations dominates the 
path integral, without being able to check that they are indeed saddle points of a full, non- 
perturbative formulation (see [7] for a recent example), or otherwise postulate a separation 
of scales between the fundamental quantum excitations of the geometry and a (larger) 
length scale characteristic for the topology changes, for example, the size of wormholes . 
We are not aware of any evidence from nonperturbative approaches that would support ei- 
ther of these assumptions. On the contrary, as already mentioned, path integrals including 
a sum over topologies tend not to exist at all. Furthermore, as is clear from studies of both 
Euclidean jH] and Lorentzian sums over geometries in two dimensions, once topology 
changes are allowed dynamically, they occur everywhere and at all scales.^ 

A new idea to tame the divergences associated with topology changes in the path inte- 

^Even in the simplest case of two-dimensional geometries, the number of possible configurations grows 
faster than exponentially with the volume of the geometry. A well-known manifestation of this problem is 
the non-Borel summability of the genus expansion in string theory. This does not necessarily mean that 
there is no underlying well-defined theory, but even in the much-studied case of two-dimensional Euclidean 
quantum gravity no physically satisfactory, unambiguous solution has been found [HI ■ 

^In the Lorentzian case cited, the standard sum over causal geometries of a fixed spacetime topology 
|1U| is extended by allowing "baby universes" which branch off the main universe, causing changes in the 
topology of spatial slices as a function of time, but not changing the spacetime topology. 
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gral was advanced in jTT] and implemented in a model of two-dimensional nonperturbative 
Lorentzian quantum gravity. The idea is to include a sum over topologies, or over some 
subclass of topologies, in the state sum, but to restrict this class further by certain geo- 
metric (as opposed to topological) constraints. These constraints involve the causal (and 
therefore Lorentzian) structure of the spacetimes and thus would have no analogue in a 
purely Euclidean formulation. In the concrete two-dimensional model considered in [TT] . 
the path integral is taken over a geometrically distinguished class of spacetimes with ar- 
bitrary numbers of "wormholes", which violate causality only relatively mildly (see also 
[12]). As a consequence, the nonperturbative path integral turns out to be well defined. 
This is an extension of the central idea of the approach of causal dynamical triangulations, 
namely, to use physically motivated causality restrictions to make the gravitational path 
integral better behaved (see [121 for a review). 

In this paper, we will present a complete analytical solution of the statistical model 
of two-dimensional Lorentzian random geometries introduced in |llj| , whose starting point 
is a regularized sum over causal triangulated geometries including a sum over topologies. 
For a given genus (i.e. number of (worm)holes in the spacetime) not all possible triangu- 
lated geometries are included in the sum, but only those which satisfy certain causality 
constraints. As shown in ^1], this makes the statistical model well defined, and an un- 
ambiguous continuum limit is obtained by taking a suitable double-scaling limit of the 
two coupling constants of the model, the gravitational or Newton's constant and the cos- 
mological constant. The double-scaling limit presented here differs from the one found 
in im IT^ . where only the partition function for a single spacetime strip was evaluated. 
We will show that when one includes the boundary lengths of the strip explicitly - as is 
necessary to obtain the full spacetime dynamics - the natural renormalization of Newton's 
constant involves the boundary "cosmological" coupling constants conjugate to the bound- 
ary lengths. Although the holes we include exist only for an infinitesimal time, and we 
do not keep track of them explicitly in the states of the Hilbert space, their integrated 
effect is manifest in the continuum Hamiltonian of the resulting gravity theory. As we will 
see, their presence leads to an effective lowering of the cosmological constant and therefore 
represents a concrete and nonperturbative implementation of an idea much discussed in 
the late eighties in the context of the ill-defined continuum path integral formulation of 
Euclidean quantum gravity (see, for example, [HI (111)- 

The remainder of the paper is structured as follows. In the next section, we briefly 
describe how a nonperturbative theory of two-dimensional Lorentzian quantum gravity can 
be obtained by the method of causal dynamical triangulation (CDT), and how a sum over 
topologies can be included. For a more detailed account of the construction of topology- 
changing spacetimes and the geometric reasoning behind the causality constraints we refer 
the reader to fU^j. The main result of Secj^lis the computation of the Laplace transform 
of the one-step propagator of the discrete model for arbitrary boundary geometries. In Sec. 
|3]we make a scaling ansatz for the coupling constants and show that just one of the choices 
for the scaling of Newton's constant leads to a new and physically sensible continuum 
theory. We calculate the corresponding quantum Hamiltonian and its spectrum, as well as 
the full propagator of the theory. Using these results, we compute several observables of 
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the continuum theory in Sec. most importantly, the expectation value of the number of 
holes and its spacetime density. In Sec. IHl we summarize our results and draw a number of 
conclusions. In Appendix A, we discuss the properties of alternative scalings for Newton's 
constant which were discarded in the main text. This also establishes a connection with a 
previous attempt [15. to generalize the original Lorentzian model without topology changes. 
In Appendix B, we calculate the spacetime density of holes from a single infinitesimal 
spacetime strip. 

2 Lorentzian sum over topologies 

Our aim is to calculate the (l+l)-dimensional gravitational path integral 



nonperturbatively by using the method of Causal Dynamical Triangulations (CDT).'^. The 
sum in denotes the inclusion in the path integral of a specific, causally preferred class 
of fluctuations of the manifold topology. The action S{g^y) consists of the usual Einstein- 
Hilbert curvature term and a cosmological constant term. Since we work in two dimensions, 
the integrated curvature term is proportional to the Euler characteristic x = 2 — 20 — 6 of 
the spacetime manifold, where g denotes the genus (i.e. the number of handles or holes) 
and h the number of boundary components. Explicitly, the action reads 



where K = l/G^ is the inverse Newton's constant and A the cosmological constant (with 
dimension of inverse length squared). 

Just like in the original CDT model [TOl (which from now on we will also refer to as 
the "pure" model, i.e. without topology changes), we will first regularize the path integral 
(H)) by a sum over piecewise flat two-dimensional spacetimes, whose flat building blocks are 
identical Minkowskian triangles, all with one space-like edge of squared length and two 
time-like edges of squared length —aa^, where a is a real positive constant. The CDT path 
integral takes the form of a sum over triangulations, with each triangulation consisting of 
a sequence of spacetime strips of height At = 1 in the time direction. A single such strip is 
a set of kn triangles pointing up and lout triangles pointing down (Figure H}. Because the 
geometry has a sliced structure, one can easily Wick-rotate it to a triangulated manifold 
of Euclidean signature by analytically continuing the parameter a to a real negative value 
|17j . For simplicity, we will set a = —l in evaluating the regularized, real and Wick-rotated 
version of the path integral (0) . 

In the pure CDT model the one-dimensional spatial slices of constant proper time t 
are usually chosen as circles, resulting in cylindrical spacetime geometries. For our present 

•^For an introduction to CDT the reader is referred to |10l 1161 ITT) . 
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Figure 1: A section of a sequence [t,t + 2] of two spacetime strips of a triangulated 
two-dimensional spacetime contributing to the regularized path integral without topology 
changes. 



purposes, we will enlarge this class of geometries by allowing the genus to be variable. We 
define the sum over topologies by performing surgery moves directly on the triangulations 
to obtain regularized versions of higher-genus manifolds [TT| I12j. They are generated by 
adding tiny wormholes that connect two regions of the same spacetime strip. Starting 
from a regular strip of topology [0, 1] x and height At = 1, one can construct a hole by 
identifying two of the strip's time-like edges and subsequently cutting open the geometry 
along this edge (Fig. Hj). By applying this procedure repeatedly (obeying certain causality 
constraints fH^]), more and more wormholes can be created. Once the regularized path 
integral has been performed, including a sum over geometries with wormholes, one takes a 
continuum limit by letting a —>■ and renormalizing the coupling constants appropriately, 
as will be described in the following sections. 

Note that our wormholes are minimally causality- and locality- violating in that they are 
located within a single proper-time step (the smallest time unit available in the discretized 
theory) and the associated baby universes which are born at time t are reglued at time 
t + 1 "without twist" fUn^- In a macroscopic interpretation one could describe them 
as wormholes which are instantaneous in the proper-time frame of an ensemble of freely 
falling observers. Note that this is invariantly defined (on Minkowski space, say) once 
an initial surface has been chosen. Such a restriction is necessary if one wants to arrive 
at a well-defined unitary evolution via a transfer matrix formalism, as we are doing. To 
include wormholes whose ends lie on different proper-time slices, one would have to invoke 
a third-quantized formulation, which would very likely result in macroscopic violations of 
causality, locality and therefore unitarity, something we are trying to avoid in the present 
model. 

One could wonder whether the effect in the continuum theory of our choice of wormholes 
is to single out a preferred frame or coordinate system. Our final result will show that this 
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Figure 2: Construction of a wormhole by identifying two time-like edges of a spacetime 
strip and cutting open the geometry along the edge. 



is not the case, at least not over and above that of the pure gravitational model without 
topology changes. The effect of the inclusion of wormholes turns out to be a rather mild 
"dressing" of the original theory without holes. We believe that the essence of our model 
lies not so much in how the wormholes are connected, because they do not themselves 
acquire a nontrivial dynamics in the continuum limit. Rather, it is important that their 
number is sufficiently large to have an effect on the underlying geometry, but on the other 
hand sufficiently controlled so as to render the model computable. 

A similar type of wormhole has played a prominent role in past attempts to devise a 
mechanism to explain the smallness of the cosmological constant in the Euclidean path 
integral formulation of /onr-dimensional quantum gravity in the continuum El • The 
wormholes considered there resemble those of our toy model in that both are non-local iden- 
tifications of the spacetime geometry of infinitesimal size. The counting of our wormholes 
is of course different since we are working in a genuinely Lorentzian setup where certain 
causality conditions have to be fulfilled. This enables us to do the sum over topologies 
completely explicitly. Whether a similar construction is possible also in higher dimensions 
is an interesting, but at this stage open question. 

3 Discrete solution: the one-step propagator 

For the (l+l)-dimensional Lorentzian gravity model including a sum over topologies, the 
partition function of a single spacetime strip of infinitesimal duration with summed-over 
boundaries was evaluated in ^1] and ^2]- In the present paper, we will extend this 
treatment by calculating the full one-step propagator, or, equivalently, the generating 
function for the partition function of a single strip with given, fixed boundary lengths. 
This opens the way for investigating the full dynamics of the model. 
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The discrete set-up described above leads to the Wick-rotated one-step propagator 

where k is the bare inverse Newton's constant and A the bare (dimensionless) cosmological 
constant, and we have omitted an overall constant coming from the Gauss-Bonnet integra- 
tion. The sum in (jH)) is to be taken over all triangulations with 1^ space-like links in the 
initial and lout space-like links in the final boundary. Note that the number of holes does 
not appear as one of the arguments of the one-step propagator since we only consider holes 
that exist within one strip. Consequently, the number of holes does not appear explicitly 
as label for the quantum states, and the Hilbert space coincides with that of the pure CDT 
model. Nevertheless, the integrated effect of the topologically non-trivial configurations 
changes the dynamics and the quantum Hamiltonian, as we shall see. 
The one-step propagator (jS)) defines a transfer matrix T by 

G \^n{lin,lout, 1) {lout\T\li'ri) , (4) 

from which we obtain the propagator for t time steps as usual by iteration, 

G X^niliriilout, t) (lout\T \lin} ■ (5) 

For simplicity we perform the sum in Q over triangulated strips with periodically identified 
boundaries in the spatial direction and one marked time-like edge. By virtue of the latter, 
G\^nilin,lout,t) satisfies the desired composition property of a propagator, 

Gx^K{lin,lout,ti + = ''^^Gx,K{lin,l,ti)Gx,K{l,lout,t2), (6) 

I 

Gx,K{lin,lout,'t + 1) = ''^^Gx,K{lin,l,^)Gx,K{l,lout,t), (7) 

I 

where the sums on the right-hand sides are performed over an intermediate constant-time 
slice of arbitrary discrete length /. 

Performing the fixed-genus part of the sum over triangulations in yields 

W2] 

-2Kg 



with N = lin + lout- To simplify calculations we will use the generating function formalism 
with 

oo 

G{x,y,g,h,l) = ^ GA,«(/m, /o«t, (9) 

where we have defined g = e~^ and h = e~'^. The quantities x and y can be seen as purely 
technical devices, or alternatively as exponentiated bare boundary cosmological constants 

x = e"^"", ?/ = e~^°'". (10) 
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Upon evaluating the sum over kn and lout one obtains the generating function of the one- 
step propagator 

1 2 

G{x,y,g,h,l) = — ■ — , ^ (11) 

1-g {x + y)l + Vl - 4^2 

with 

h , , 

1 ^ - (12) 

Note that in order to arrive at the final result (fTT|) . we have performed an explicit sum 
over all topologies! The fact that this infinite sum converges for appropriate values of the 
bare couplings has to do with the causality constraints imposed on the model, which were 
geometrically motivated in JTj, and which effectively reduce the number of geometries in 
the genus expansion. 

In fjllj) one recognizes the generating function CsX{u^^ for the Catalan numbers, 

For h = one has Cat{u'^) = 1 and expression (fTTj) reduces to the one-step propagator 
without topology changes, 

G{x, y,g,h = 0, 1) = ^— -. (14) 

I- g [x + y) 

Furthermore, one recovers the one-step partition function with summed-over boundaries 
of (HUE] by setting x = y = l, 

Zig,h,l) = Y^ , ^ ^ • (15) 

1-2^1 4(^)2 

4 Taking the continuum limit 

Taking the continuum limit in the case without topology changes is fairly straightforward 
[TU] . The joint region of convergence of ()14j) is given by 

\x\ < 1, \y\ < 1, 1^1 < ^. (16) 
One then tunes the couplings to their critical values according to the scaling relations 

g = ^{l-a'A) + 0{a'), (17) 
x = l-aX + 0{a^), y = l-aY + 0{a'^). (18) 

Up to additive renormalizations, x, y and A scale canonically, with corresponding renor- 
malized couplings X, Y and A. In the case with topology change we have to introduce an 
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additional scaling relation for h. Since Newton's constant is dimensionless in two dimen- 
sions, there is no preferred canonical scaling for h. We make the multiplicative ansatz^ 



h = l^hren{adf, (19) 

v2 

where hren depends on the renormalized Newton's constant Gat according to 

Ken = e-^^"'^. (20) 

In order to compensate the powers of the cut-off a in ()19|) . d must have dimensions of 
inverse length. The most natural ansatz in terms of the dimensionful quantities available 
is 

d= (VA"(x + r)^-"). (21) 

The constants (3 and a in relations ()19p and ()2H) must be chosen such as to obtain a 
physically sensible continuum theory. By this we mean that the one-step propagator should 
yield the Dirac delta-function to lowest order in a, and that the Hamiltonian should be 
bounded below and not depend on higher-order terms in (fTTj) . (fT^ . in a way that would 
introduce a dependence on new couplings without an obvious physical interpretation. 

To calculate the Hamiltonian operator H we use the analogue of the composition law 
(P) for the Laplace transform of the one-step propagator [TU] . 

/dz 
-—G{x,z-'-l)G{z,y,t). (22) 

In a similar manner we can write the time evolution of the wave function as 

V^(x, t + l) = j ^Gix, z-'; l)^{z, t). (23) 

When inserting the scaling relations ()17|) . (fTSj) and t = — into this equation it is convenient 
to treat separately the first factor in the one-step propagator (fTT|) . which is nothing but 
the one-step propagator without topology changes (fT^ . and the second factor, the Catalan 
generating function (fTSj) . Expanding both sides of (PHj) to order a gives 

{l-aH, 0(a^)) m) - £ i { + « Cat(,.^)} (24) 

where we have used 

^(X,T + a) = e-"^V'(^,T), (25) 

with ip{X) = ilj{x = l — aX). Note that the first term on the right-hand side of ()24j) . is 
the Laplace-transformed delta-function. The interesting new behaviour of the Hamiltonian 



"^Here the factor is chosen to give a proper parametrization of the number of holes in terms of 
Newton's constant (see Section 13)). 
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is contained in the expansion of the Catalan generating function. Combining ()13p and ()12|). 
and inserting the scahngs (fTTj) . (fTSj) . yields 



Cat{u') = l + j^^a^'-^ + h.o., (26) 

where h.o. refers to terms of higher order in a. In order to preserve the delta-function 
and have a non-vanishing contribution to the Hamiltonian one is thus naturally led to 
P = 3/2. For suitable choices of a it is also possible to obtain the delta-function by setting 
(3=1, but the resulting Hamiltonians turn out to be unphysical or at least do not have an 
interpretation as gravitational models with wormholes, as we will discuss in Appendix A.^ 
For (3 = 3/2 the right-hand side of (j24j) becomes 



too 



2m\Z-X \{X-Zy (X - zf 



We observe that for a ^ the last term in (j27p does not contribute to the Hamiltonian. 
Performing the integration for a > and discarding the possibility of fractional poles the 
Hamiltonian reads 

HiX, ^) = X2 A + X - 2A A + 2A^/i,\, ^~/^'°5r\ , a = -, -, 1, ... (28) 
^ 'dx' dX dX r(3«) ' 3'3' ' v ; 

For all a's, these Hamiltonians do not depend on higher-order terms in the scaling of the 
coupling constants. One can check this by explicitly introducing a term quadratic in a 
(which can potentially contribute to H) in the scaling relations (jl8p . namely, 

X = 1-aX + ^-fa^X^ + 0{a^), 

y = i-aY + ^ja^Y^ + 0{a^), (29) 

and noticing that (j^Hj) does not depend on 7. After making an inverse Laplace transforma- 
tion ■?/'(L) = dXe'^ ^ip{X) to obtain a wave function in the "position" representation 
(where it depends on the spatial length L of the universe), and introducing m = 3a — 1 
the Hamiltonian reads 

d 8"^ 8 2A'^/?2 

*<^'al)--^aIJ-aI + ^^^-T(;^^" '-".i.^.-- (30) 

Since H is unbounded below for m ^ 2, we are left with m = and m = 1 as possible choices 
for the scaling. However, setting m = merely has the effect of adding a constant term to 



^One might also consider scalings of the form h cihrenio.d)+C2hreniO'd)^^^ , but they can be discarded 
by arguments similar to those of Appendix A. 
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the Hamiltonian, leading to a trivial phase factor for the wave function. We conclude that 
the only new and potentially interesting model corresponds to the scaling with m = 1 and 



/^2 = l/.LA(X + r)a^ (31) 



with the Hamiltonian given by 



^) = - ^ + (1 - h'en) 2 A L. (32) 

Note that for all values Gat > of the renormalized Newton's constant (|19|) the Hamiltonian 
is bounded from below and therefore well defined. It is self-adjoint with respect to the 
natural measure dfi{L) =dL and has a discrete spectrum, with eigenfunctions 



ML) = A„e-V2A{i-ejii:^(2^2A(l - hlJ ), n = 0, 1, 2, ... , (33) 
where L„ denotes the n'th Laguerre polynomial. Choosing the normalization constants as 



X= v/8A(l-/i,U, (34) 

the eigenvectors {tpn{L), n = 0, 1,2, ...} form an orthonormal basis, and the corresponding 
eigenvalues are given by 

E„ = v/2A(l - J (2n + l), n = 0,l,2,.... (35) 

Having obtained the eigenvalues one can easily calculate the Euclidean partition function 
for finite time T (with time periodically identified) 

°o -V2A(l-/x2,„)T 

For completeness we also compute the finite-time propagator 

oo 

Ga,gALuL2,T) = (L2|e-^^|Li) = e-^^"C(L2)^„(Li). (37) 

n=0 

Inserting ()33p into ()37|) and using known relations for summing over Laguerre polynomials 
[ini yields 



Ga,g. (Li, L,, T)-u -^^^ Jo j 



(3J 



where we have used the shorthand notation uj = ^J2K{1 — h'^.^n)- As expected, for hyen 
the results reduce to those of the pure two-dimensional CDT model. 
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5 Observables 

Due to the low dimensionality of our quantum-gravitational model, it has only a few 
observables which characterize its physical properties. Given the eigenfunctions (jH!^ of 
the Hamiltonian (jH^ one can readily calculate the average spatial extension (L) of the 
universe and all higher moments 

POO 

(L-)„= / dLL"^\ML)\'. (39) 
Using integral relations for the Laguerre polynomials pSJ one obtains^ 

m 

,rm^ ( 1 V r(n-m)r(m+i) 



^8A(l-/i2^j; r(n + l)r(-m) 
X 3F2(— n, 1 + m, 1 + m; 1, 1 + m — n; 1), (40) 
where 3-^2(01, 02, as; 61, 62; z) is the generalized hypergeometric function defined by 



00 



sF,{a,, a2, as; 61, h; z) = ^ (bMhUl " ^''^ 

Observe that the moments scale as {L'"^)^ ~ which indicates that the effective Haus- 
dorff dimension is given by dn = 2, just like in the pure CDT model 

In addition to these well-known geometric observables, the system possesses a new type 
of "topological" observable which involves the number of holes A^'g, as already anticipated 

^Note that the poles of r(— m) cancel with those of the hypergeometric function. 
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in fn\ \1'2\. As spelled out there, the presence of holes in the quantum geometry and their 
density can be determined from light scattering. An interesting quantity to calculate is 
the average number of holes in a piece of spacetime of duration T, with initial and final 
spatial boundaries identified. Because of the simple dependence of the action on the genus 
this is easily computed by taking the derivative of the partition function Zt with respect 
to the corresponding coupling, namely, 

(ATg) = J_^M^_ (42) 
Zt 2 dhyen 

Upon inserting (jH^ this yields 



(v/2A(T 



coth U/2Afl-/i2 )T 



ren/ 



<^°>^^"-" V2A(l-/4„) ^ '''' 

In an analogous manner we can also calculate the average spacetime volume 

leading to 

{V) = T ^^^Zj'-^ coth [Vmi-hUT) . (45) 



/2A 

Dividing ()43|) by ()45|) we find that the spacetime density n of holes is constant, 

{V) I -Ken ^ ^ 

The density of holes in terms of the renormalized Newton's constant is given by 

n=^ A. (47) 

e'^iv — 1 

The behaviour of n in terms of the renormalized Newton's constant is shown in Fig. 01 The 
density of holes vanishes as Gat —> and the model reduces to the case without topology 
change. - An alternative calculation of the density of holes from an infinitesimal strip, 
which leads to the same result, is presented in Appendix B. 

We can now rewrite and interpret the Hamiltonian ()32j) in terms of physical quantities, 
namely, the cosmological scale A and the density of holes in units of A, i.e. rj = j^, resulting 
in 

d d"^ d 1 

One sees explicitly that the topology fluctuations affect the dynamics since the effective 
potential depends on r], as illustrated by Fig. HI 

It should be clear from expressions ()48|) and (jTfj) that the model has two scales instead 
of the single one of the pure CDT model. As in the latter, the cosmological constant defines 
the global length scale of the two-dimensional "universe" through (L) ~ The new scale 
in the model with topology change is the relative scale t] between the cosmological and 
topological fluctuations, which is parametrized by Newton's constant Gn- 
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Figure 4: The coefficient of the effective potential, f{r]) = 1/(1 + rj), as function of the 
density of holes in units of A, rj = j. 



6 Conclusions 

In this paper, we have presented the complete analytic solution of a previously proposed 
model JI] of two-dimensional Lorentzian quantum gravity including a sum over topologies. 
The presence of causality constraints imposed on the path-integral histories - physically 
motivated in fT| |^ - enabled us to derive a new class of continuum theories by taking an 
unambiguously defined double-scaling limit of a statistical model of simplicially regularized 
spacetimes. After computing the Laplace transform of the exact one-step propagator of the 
discrete model, we investigated a two-parameter family, defined by (fT^ and (|7T|). of possible 
scalings for the gravitational (or Newton's) coupling, from which physical considerations 
singled out a unique one. For this case, we computed the quantum Hamiltonian, its 
spectrum and eigenfunctions, as well as the partition function and propagator. Using 
these continuum results, we then calculated a variety of physical observables, including the 
average spacetime density of holes and the expectation values of the spatial volume and 
all its moments. 

This should be contrasted with the previous treatment in fUE], in which only the one- 
step partition function with summed-over boundaries was evaluated. Because of the lack of 
boundary information, no explicit Hamiltonian was obtained there. Moreover, it turns out 
that the scaling of the couplings which in the current work led to the essentially unique 
Hamiltonian (jH^ could not have been obtained or even guessed in the previous work. 
This is simply a consequence of the fact that the dimensionful renormalized boundary 
cosmological constants make an explicit appearance in the scaling relation (|31|) for h, 
and thus for Newton's constant. We conclude that - unlike in the case of the original 
Lorentzian model - for two-dimensional causal quantum gravity with topology changes 
one cannot obtain the correct scalings for the "bulk" coupling constants from the one-step 
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partition function with boundaries summed over (which is easier to compute than the full 
one-step propagator). 

In contrast with what was extrapolated from the single-strip model in ^2], the total 
number of holes in a finite patch of spacetime turns out to be a finite quantity determined by 
the cosmological and Newton's constants. Note that this finiteness result has been obtained 
dynamically and without invoking any fundamental discreteness. Since the density of holes 
is finite and every hole in the model is infinitesimal, this implies - and is confirmed by 
explicit calculation - that the expectation value of the number of holes in a general spatial 
slice of constant time is also infinitesimal. The fact that physically sensible observables 
are obtained in this toy model reiterates the earlier conclusion ^1] that causality-inspired 
methods can be a useful tool in constructing gravitational path integrals which include a 
sum over topologies. 

From the effective potential displayed in Fig. |3 one observes that the presence of worm- 
holes in our model leads to a decrease of the "effective" cosmological constant /(?7)A. In 
Coleman's mechanism for driving the cosmological constant A to zero [HlEl; an additional 
sum over different baby universes is performed in the path integral, which leads to a distri- 
bution of the cosmological constant that is peaked near zero. We do not consider such an 
additional sum over baby universes, but instead have an explicit expression for the effective 
potential which shows that an increase in the number of wormholes is accompanied by a 
decrease of the "effective" cosmological constant. A first step in establishing whether an 
analogue of our suppression mechanism also exists in higher dimensions would be to try 
and understand whether one can identify a class of causally preferred topology changes 
which still leaves the sum over geometries exponentially bounded. We will return to this 
issue in a future publication. 
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A Scalings with /3 = 1 

In this appendix we discuss the scalings with (3=1 which we discarded as unphysical in 
Sec. m above. We proceed as before by inserting the scaling relations (fTTj) and into 
the composition law (f^ . Instead of using /3 = | we set /5 = 1, leading to the scaling 

h='^hrenaVA'^{X + YY-'', (49) 
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where the normahzation factor on the right-hand side has been chosen for later convenience. 
Up to first order in a one obtains 

/too ly 
— {A{X, Z) + B{X, Z)a + 0{a^)]iP{Z), (50) 
-ioo ^'^'^ 

where the leading-order contribution is given by 



■^'^•^>= (Z-.Y)(l + C(A-.Z)) <^^' 



with 



c(x, z) = ^i^h^Jx^zy^- (52) 

For the Laplace transform of A(X, Z) to yield a delta-function, the scaling should be chosen 
such that a ^ 0. Considering now the terms of first order in a, 

BiX Z) = Ke^jX + Z - AZ^)A- 



(X - z)i+2"C(x, z) (1 + c{x, z)y 

XZ -2A + j{X -Zf 

{x-zyc{x, z){i + c{x, z))' 



(53) 



one finds that for a ^ — 1 the continuum limit is independent of any "hole contribution" 
(i.e. terms depending on hren) and therefore leads to the usual Lorentzian model. This 
becomes clear when one expands the last term of in (X — Z), resulting in 



(X-Z)2C(1 + C) 1{X-ZY\ 4 

For a ^ — 1 the term depending on hren does not have a pole and therefore does not 
contribute to the Hamiltonian. Since we are only interested in non-fractional poles, this 
leaves as possible a-values only a = and a = —\. 



A.l The case = 1, a = 

For a = the Hamiltonian retains a 7-dependence contained in the first line of ()53p. Since 
there is no immediate physical interpretation of 7 in our model, it seems natural to choose 
7 = 0, although strictly speaking this does not resolve the problem of explaining the 7- 
dependence of the continuum limit. Setting this question aside, one may simply look at 
the resulting model as an interesting integrable model in its own right. In order to obtain 
a delta-function to leading order, one still needs to normalize the transfer matrix by a 
constant factor 2/(1 + s), with s := a/1 — /i^g^. After setting 7 = and performing an 
inverse Laplace transformation, the Hamiltonian reads 

(9 1 / 92 d \ 
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It is self-adjoint with respect to the measure djj,{L) = ^dL. Further setting L = ^ one 
encounters the one-dimensional Calogero Hamiltonian 

- . d , 1 92 1 1 A 

with u = and A = l — 4(1 — s)^, which implies that the model covers the parameter 
range —3 ^ A ^ 1. The maximal range for which the Calogero Hamiltonian is self-adjoint 
is —00 < A ^ 1. The usual Lorentzian model without holes corresponds to ^4 = 1. The 
Hamiltonian ()56|) has already appeared in a causal dynamically triangulated model where 
the two-dimensional geometries were decorated with a certain type of "outgrowth" or small 
"baby universes" ^H]- This model covered the parameter range ^ A ^ 1. 
The eigenvectors of the Hamiltonian are given by 



^n{L) = AnC-^^ iFi{-n, s, 2V2AL), dfi{L) = L'-^dL, (57) 

where iFi{—n,a,b) is the Kummer confluent hypergeometric function. The eigenvectors 
form an orthonormal basis with the normalization factors 



and the corresponding eigenvalues 



v2A 

En = - (2n + s), n = 0,1,2,.... (59) 

s 

One sees explicitly that the case s = 1 or, equivalently, A = 1 corresponds to the pure 
two-dimensional CDT model. 

A. 2 The case (3 = I, a = -\ 

For a = — I the result does not depend on 7 and therefore on the detailed manner in which 
we approach the critical point. However, the Hamiltonian 

cannot be made self-adjoint with respect to any measure dfi{L) because the boundary part 
of the partial integration always gives a nonvanishing contribution. We therefore discard 
this possibility. 
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B The density of holes of an infinitesimal strip 



In this appendix we give an alternative derivation of the spacetime density n of holes and 
explicitly show that the number of holes in a spacetime strip of infinitesimal time duration 
a is also infinitesimal. The operator in the L-representation of the number of holes per 
infinitesimal strip with fixed initial boundary L can be calculated by 

N,,a-.o = T-'^-—, (61) 

where T is the transfer matrix defined in Using T = 1 — aH + 0{a?) and evaluating 
flUTj) to leading order in a gives 

N,,a^o = -a^^ + 0{a') = 2Ahl^La + 0{a'). (62) 

Similarly, the volume operator of the same infinitesimal spacetime strip in the L- representation 
is given by 

rll^ f)f-f 

Va-,0 = -T-'^ = ^ + = 2 (1 - hln)La + 0{a'). (63) 

Although both expressions (jU^ and (jUHj) vanish in the limit as a — *■ (and therefore the 
number of holes and the strip volume are both "infinitesimal"), their quotient evaluates to 
a finite number independent of L, namely, 

n = ^^ = -^^A. (64) 

This is the exactly the same result for the spacetime density n of holes as we obtained 
earlier from the continuum partition function (j46|) . 
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